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Abstract 



The search for particle electric dipole moments (edm) represents a 

most promising way to search for physics beyond the standard model. 

A number of groups are planning a new generation of experiments using 

stored gases of various kinds. In order to achieve the target sensitivities 

it will be necessary to deal with the systematic error resulting from the 

. — > 

interaction of the well-known v x E field with magnetic field gradients 
(often referred to as the geometric phase effect (Commins, ED; Am. J. 
Phys. 59, 1077 (1991), Pendlebury, JM et al; Phys. Rev. A70, 032102 
(2004)). This interaction produces a frequency shift linear in the electric 
field, mimicking an edm. In this work we introduce an analytic form 
for the velocity auto-correlation function which determines the velocity- 
position correlation function which in turn determines the behavior of 
the frequency shift (Lamoreaux, SK and Golub, R; Phys. Rev A71, 
032104 (2005)) and show how it depends on the operating conditions of 
the experiment. We also discuss some additional issues. 
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1 Introduction 

The proposition that the search for particle electric dipole moments (edm) rep- 
resents a reasonable method to look for physics beyond the standard model PP 
is inspiring many groups to search for edm's in a variety of systems. (See .2 for 
a recent review). Experiments on several systems including the neutron [3J, and 
several species of confined gases 4 including Radium 5 , Radon |S] and Xenon 
are in various stages of preparation. These experiments are all hoping to 
reach sensitivities in the range of 10 -27 — 10~ 28 e — cm. Sensitivity in this range 
has already been achieved in the case of Hg The experiments proposed 
represent a broad range of operating conditions, from room temperature gases 
with buffer gas to laser cooled atoms in a MOT. 

In order to achieve the target sensitivities it will be necessary to deal with 
the systematic error resulting from the interaction of the well-known 1/ x E field 
with magnetic field gradients. Often referred to as the geometric phase effect 
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IS], EH this interaction produces a frequency shift linear in the electric field, 
mimicking an edm. This systematic effect is highly dependent on the operating 
conditions of the experiment. While experiments in small vessels and with high 
pressure buffer gas are expected to be relatively insensitive to the systematic 
effect, each of the proposed experiments will have to be analyzed in detail to 
judge its sensitivity to the effect and to find methods of dealing with it. 

In this work we introduce an analytic form of the correlation function which 
determines the behavior of the frequency shift JT] and show in detail how it 
depends on the operating conditions of the experiment. For clarity we special- 
ize the discussion to the Los Alamos proposal for a neutron edm search using 
Ultra-cold neutrons (UCN) and He 3 atoms diffusing in supcrnuid He 4 as a 
co-magnetometer, ^2] but the generalization to other cases is straightforward. 

First analyzed by Commins |5] in the context of a beam experiment, the 
frequency shift has been discussed in some detail by Pendlebury et al ^Oj in 
connection with experiments involving stored particle gases. Additional discus- 
sion and calculations have been given by 

Our present understanding of the effect can best be summarized by figure 
1, which appeared as figure 3 in This is a plot of the normalized (linear in 
E) frequency shift, 5u>, vs. normalized Larmor frequency loq, for various values 
of collision mean free path, A, and wall specularity, calculated for particles 
moving with fixed velocity, v, in a cylindrical measurement cell of radius R. 
The horizontal scale is fixed by the frequency of motion around the cell (v/R). 
In general the shift for UCN will be given by a value of ui / (v/R) > 4 while 
co- magnetometer atoms are characterized by lu / (v/R) <C I. 

These results have been obtained by numerical simulation of the position- 
velocity correlation function and taking the Fourier transform. According to 
[fTT]. Eq. 26] the frequency shift is given by 



where a — ^dB z /dz, b = jE/c, 7 is the gyromagnetic ratio, and R{t) is the 
position-velocity correlation function for motion in the plane perpendicular to 
the z axis, defined in [|^, Eq.. 27]: 



From the experimental point of view it is very appealing to try to make use of 
the zero crossing, apparent in figure I, to reduce the effect. 

In this note we present an analytic form for the velocity correlation function 
from which R (t) can be determined and compare it to the results obtained 
previously by numerical simulations. In the collision- free case we obtain the 
result obtained in (|TJ]|, Eq.. 78) for a single trajectory. Using the analytic 
function for the case with gas collisions, we average over a Maxwell velocity 
distribution and calculate the temperature dependence of the frequency shift 
for 3 He diffusing in superfluid 4 He. We also propose a method for measuring 



T 




CD 







R (r) = (rl(t) • vj(t - r) - r±(t - r) • v±(t)) . 



(2) 
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the spectrum of the correlation function, i. e. the frequency dependence of the 
shift, directly. 



2 Analytical form for the correlation function 

R(r) 

According to ^I] the correlation function R(t) is determined by the velocity 
autocorrelation function, 



namely, 



Notice, that 



ij(t) = (v x (t)-v x (0)), (3) 

T 

R( T ) = 2 / 4>(t)dt. (4) 



o 



R{t) — * 0, when r — * oo. (5) 
Thus we start by considering ij)(t). 

2.1 Specular wall collisions, no gas collisions 
2.1.1 Specification of trajectories 

We consider particles moving ballistically in a cylindrical storage cell with a 
fixed velocity v. As shown in JOIjEB the frequency shift depends only on the 
motion in the x, y plane. Referring to figure 2, the trajectory sweeps out an 
angle 

a = arccos (6) 
with respect to the center in a time 

2i?sina . , 

Tw = 7 

V 

where t w is the time between wall collisions. For specular reflections the angle 
a, characterizing the trajectory and the velocity are unchanged by reflection. 
The average angular velocity for a single trajectory is then 

u(a) = — . (8) 

Let F(r) dr be the probability that the trajectory of a particle has a distance 
of closest approach to the center in the interval (r, r + dr). According figure 2, 



2VR 2 - r 2 dr A^R 2 - r 2 
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We will use the distribution P(a)da of the trajectories over the angle a, where 

tt/2 

rlr A Qin z n 

P(a)|da| = F(r)\dr\ P(a) = F(r) 



dr 
da 



J P{a)da = 1, 
o 

(10) 



a result obtained in |1U) . 
2.1.2 The velocity autocorrelation function 

We now calculate the velocity autocorrelation function for the particles moving 
along the trajectories with given a (or the pericenter r = i?cosa), 

i;(a,t) = (vx(t)-ifx(0)), ^(a,0) = «i. (11) 

Here the averaging goes only over the initial position of the particles. Let the 
velocity autocorrelation function be denoted by 

_ v_i_{t) -i/j.(0) . . 

f{x, a, t) = - . (12) 

for a particle on a trajectory characterized by a, starting at the position x, 
measured from the end of a chord, at t = 0. Thus, 

2i? sin a 
V 2 [ 

^( a '*) = 7T5~- / f(x,a,t)dx. (13) 

zKsina J 
o 

As the speed \v \ = v is not changed by the collisions the velocity correlation 
function is proportional to cos#(t) where 9(t) is the angle between v±(t) and 
v± (0) . The starting position x determines the exact times at which the collisions 
occur, the time between collisions being given by (JJJ in all cases. 

As a result, for the time I ■ t w < t < (I + 1) • t w (I = 0, 1, 2 . . .) the function 
/ is given by 

/(x,a,t) = | 008 J? l (i + 1)a) » °^ X , <v{ l~ l " rw l , w ■ (M) 
J v ' ' ; \ cos(2Za), v(t- I ■ t w ) < x < 2Rsma, v ' 

Then performing the averaging (|13|1 the autocorrelation function takes the form: 

i>{a,t)=v 2 {A l + B l — ), (15) 

where 

Ai = (I + 1) cos(2Za) - I cos(2(/ + l)a), (16) 
Bi = cos(2(Z + l)a) - cos(2Za). (17) 
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2.1.3 Spectrum of the velocity correlation function 

The autocorrelation function can be written as a Fourier integral, valid for both 
closed, periodic orbits and general open ones: 



ip(a,t) = / *(a, u>) cos(iot)du), 



where ty(a,uj) — ^*(a,u) and 4 , (a,w) = ^(a, — 10), so that 



ip(a, t) cos(u>t)dt. 



Straightforward calculation (see Appendix A) gives: 

2a — 2-kvtC 



(18) 



(19) 



* (a, a;) 



r> 2 • 2 

2v sin a 



E 



+ 6[u- 



2a + 2nm 



2v 2 sin 2 a 



E [ s +a(w-w-(a))] 



(20) 



where 



2 (7rm ± a) 



(21) 



Notice that these frequencies are the resonant frequencies found previously in 
the behavior of the the frequency shift as function of the angle a for particles 
moving inside a cylindrical cell without any damping - see figure 7 and equation 
(78) in [TJ]|. (Below we rederive this equation in the frame of our approach and 
show how it can be generalized to account for damping). 

2.1.4 Solution for the frequency shift in the absence of gas collisions 

Now the frequency shift Auj (a) see equation (40) ref. 
* (a,ui) 



f 

-Auj (a) — ab 



dio 



2abv 2 sin 2 a ^ / 



5 (lu + 2 % 2 J m ) + s(w- 2a+2 J m ) 



(t w u) (UJ 2 - LU 2 ) 



■dio 



2abv 2 sin 2 a 



E 



2 

T„„ X 



m— — oo 



(2a - 27rrn) 2 (t2 w 2 _ ( 2a - 27rm) 2 ) (2a + 27rm) 2 (t^u 2 , - (2a + 27rm) 2 ) 

(22) 
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The sum is over all to, so terms with ±m are included twice. Using t w = 
2Rsiria/v and writing lo t w — 2u> Rsuia/v — 2oj' sin a = 25 with uj' = uj R/v 
being the dimensionless frequency we then find: 

00 1 1 

-Auj (a) = R 2 ab sin 4 a V ^7 C ( 23 ) 

^ (a + 7rm) fa _ ( Q + ^ TO ) 2 j 



m— — 00 



^i? 2 a& sin 4 a x. 



2<5 D 



m— — 00 

1 1 \ 1 



((<5 Q - a) - 7rm) ((S + a) + nm) / (a + nm) 



(24) 



00 ^ ^ 
-Auj (a) = R 2 ab sin 2 a > . 7 ^, (25) 

™ (a + - ) 

We now evaluate (|24|) using equation (|72(1 derived in Appendix B. 

+00 



E 



f f cot 9 — cot </? 



n— — 00 



(7m - ip) 2 (-nn - 9) (ip-9)sin 2 ip (ip - 9) 2 

Rewriting JUJ): 



-Auj (a) = —R 2 ab sin 4 aA 
2o D 



with 

a- E 



m— — 00 



((tf - a) - 7rm) ' ((5 + a) + 7rm)y ( a + 7r m ) 2 
1 cot (8 — a) + cot a 1 cot (<5 Q + a) — cot a 



5 sin 2 O! S 2 <5 sin 2 a S 2 

2 cot (S — a) + cot (6 Q + a) 

S sin 2 a S 2 

2 ( sin 2 a sin 2S Q 



S sin a \ 2S Q sin ((5 D - a) sin (<5 C + a) 



we find 



, n / « \ 2 , /, sin 2 a sin 2<5 G \ ,„„. 

-Aw (a = — aft 1 + 9A . T , ; . ° - > (26) 

\w / V 2d sin (d — a) sin (o c + a)/ 

This formula was originally derived in |10j . (equ. 78) by direct solution of 
the classical Bloch equations and the result shows the equivalence of the two 
methods. 
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2.2 Influence of Gas collisions 

2.2.1 The collision- free velocity correlation function as a sum of har- 
monic oscillators 

Substituting i|2U[l into l|18|) , we obtain for the velocity autocorrelation function: 



ip(a,t) 



2v 2 sin 2 a 



E 



n— — oo 



mj+ (a)t cosw„ (a)t 



(tot (a)) (w„ (a))' 



2 sin a 



E 



(ujn (a)) 2 (w B (a)) 2 . 



(27) 
(28) 



i.e. a sum of oscillating terms (tp^ (a,t) — v 2 cos lu^ (a) i) with different frequencies. 
Each term obeys an equation: 



dt 2 



+ (uj^(a)) 2 ^(a,t) = 0. 



(29) 



Notice that the fundamental frequency (a) coincides with (jHJ). The cor- 
responding term obviously dominates in the decomposition l|28|) . 

2.2.2 Velocity correlation function in the limit of short gas-collision 
times 

We consider a particle that moves among scattering centers. If r c is the average 
time between collisions the velocity autocorrelation function will have the form 

ip(t) = (v(t)v(0)} =v 2 e~t. (30) 
In the other words, ip(t) obeys the equation: 

« + -m = o. (3D 

at t c 

2.2.3 Combined influence of gas and specular wall collisions 

The velocity autocorrelation function in the absence of collisions is given by the 
sum of the motions of a group of harmonic oscillators (|28|l . In the presence 
of gas collisions the individual oscillators tp^(a,t) will obey the equation for a 
damped harmonic oscillator which is the combination of (|31|) and l|29|) . i.e. 



dt 2 



dt 



+ u 2 Mt) = o, 



(32) 



with the initial condition: 



^t(a 1 0)=v 2 . 



(33) 
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The boundary conditions (@J and © satisfied by ip(t) mean that 

T 

J ip n (&,t)dt — > 0, when r — > oo. (34) 



Generally, the equation for a damped harmonic oscillator (|32[1 has the general 
solution: 

^{t) = de-^ 1 + c 2 e~^ t , (35) 

where 



^ = ^ c +^-" 2 > (36) 
Taking into account boundary conditions (|3^|l and l|34|). we get: 

= f e"*' - ^e-^V (37) 

m - m \ 771 ) 

The correlation function Q takes the form (Eq. 36, 

i?„(a,r) = 2 / ifj n (a,t)dt = -^(l-e-^-^ T ) e~^\ 
o 

The last expression can be rewritten in the form: 

S{LO,T c ) \ J 

with 

A = vt c , s{u, t c ) = v/l-4w 2 r c 2 . (39) 

We then have for a single oscillator (see equation 

t 

Su) n = — lim / R n (r) cos(o;ot)c?t = abS n ,d (u) ) (40) 

Z t— >oo J 



C / N 2 (w2-W») 

<Vd K) = -v 



((w*-w*) a + ^rg) 

where we again introduced wj, = u> R/v and r D = R/X c {v) with A c («) the 
velocity dependent mean free path. The frequency shift will then be a sum of 
such terms for each of the oscillators as in (|25|l . 
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Comparing to the collision-free case in 12511 we see that in the presence of 
damping the frequency shift will be given by 



-Au> (a) = i? 2 a& sin 2 a ^ 

'~ — ' (a + irm) 



m— — oo 



f a _ (a+TTm) 2 \ 
\ o sin^ a J 



( /2 _ (a+irm) 2 \ ' 
\ o sin a a J 



(42) 

that is we go from the collision free case to the case of gas collisions by replacing 

fa («') = ' 



( u> 2 - {a+ T )2 ) 

V o sin^ a J 

in (|2"5*|l by the square bracket in (|4*2"|l or by replacing /„ (u)') by f a (w'-^/l + i'^) 
and taking the real part. Since we have evaluated the summation Ij25(l we obtain 
the frequency shift by making the equivalent transformation to (|26|l 



-Aw(a) = iFabsu? aRe iF P (a,5 = 5 o\j l + ^r)\ ( 43 ) 
where 

/ sin 2 a sin 25 \ 1 

Fp(a ' ^={ 1+ 26 sin (5 -a) sin (5 + a) ) P 

(remember 6 — lu t w /2). For a fixed velocity we average over a, according to 

t/2 



Aw= / daP (a) Aw (a) 
Jo 



(44) 



The results are shown in fig. 3 in comparison with the results of the numerical 
simulations obtained in [FT]. 

We see the agreement is quite good, within the uncertainties of the numerical 
simulations. The agreement in the region of the zero-crossings is excellent. 



3 Frequency shift averaged over velocity distri- 
bution, temperature dependence 

In the neutron edm experiment proposed by the EDM collaboration |12j a dilute 
solution of He 3 dissolved in superfluid He 4 will be used as a co-magnetometer 
to monitor magnetic field fluctuations. As such the He 3 will see essentially the 
same magnetic and electric fields as the neutrons and will be subject to the 
linear E field systematic under discussion. 

Using the analytical form of the correlation function we can average the E 
field proportional frequency shift over the Maxwell velocity distribution for a 
gas in thermal equilibrium. We take the realistic case of the mean free path for 
collisions proportional to velocity, (collision time r c , independent of velocity). 
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corresponding to a cross section ~ 1/v. This applies to He 3 in superfluid 
He 4 . Since the velocity of the He 3 is much less than the phonon velocity (2.2 x 
10 4 cm/ sec) the collision rate of phonons with the He 3 will be independent of 
the He 3 velocity. Thus in a time r c a He 3 with velocity v, will move a distance 
X V (T) = vt c (T). We will obtain the collison time, t c (T), from the measured 
values of the diffusion constant for He 3 in superfluid He 4 [T5| : 

£>(T) = ^ cm 2 / sec 

Then r c (T) = 3D (T) / (v 2 ) T where (w 2 ) T is the mean square velocity in a 
volume of gas. r a is defined as 

R R R r a 



X V (T) vt c (T) y0(T)r c (T) y 

with y = v/(3(T) and (3(T) is the most probable velocity in a volume. R the 
radius of the cylindircal vessel is taken as R = 25cm in the numerical calcula- 
tions. Both (v 2 ) T and 0(T) are calculated using the effective mass of He 3 in 
the superfluid: — 7.2 amu. 

For a single velocity the frequency shift will be given by (|42|) or (|43|l averaged 
over a. Using Eq. (fTTTf) 

/■7T/2 



Alj= daP (a) Alo (a) 

Jo 



R 2 ab— / da sin 4 a 



7T 



o m (a + nmY 



, ,/2 (a+irmy 

LO „ : T 

° sin a 



/9 (a+TTTa) 2 \ . »o 9 



(45) 



We now replace w' a = uj q R/v = L0 o R/y(3(T) = Lo*/y where lo* = L0 o R/f3{T). 
Then 



4 r /2 . i 

-Aw(y) = R 2 ab— / da sin 4 a \^ 

n Jo m (a + nm) 



(a+irm) 2 2) ,.2 



5 ' y y 

sin z a ^ i a 



2 

("+7T"t) 2 „.2 1 \ _|_,.,*2~2 



where r a = r y — R/X c . Averaging over the two dimensional velocity distribu- 
tion (tii) we obtain for the normalized frequency shift 

* K ' T) = ablv J v e ~ y2 Auj (y) d y ( 46 ) 

The results are plotted in fig. 4 which shows the frequency shift as a function 
of reduced frequency lo* for various temperatures. 

In figure 5 we show an expanded plot of the normalized frequency shift in 
the region near the zero crossings as a function of temperature for fixed lo* (T). 

It is evident that the collisional damping can lead to large reductions in the 
effect for He 3 . 
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4 Non-specular wall collisions 



In this section we give a brief description of how non-specular wall collisions can 
be included in the calculation. Detailed study of this problem will be left for a 
future work. 

In the preceding sections we have shown that the velocity autocorrelation 
function can be regarded as the result of the sum of harmonic oscillators of 
different frequencies. 

Considering one such oscillator during one traversal of the cell the oscillator 
will undergo a phase change 

<j) = ojT w = 2a. (47) 

A non-specular reflection from the wall would result in a change in the incident 
angle for the next collision, x, by a random amount A% and hence a change in 
the accumulated oscillator phase by 

A0 = 2A X , (48) 

because x = § — a - 

Since the changes A(jj are random the phase (f> will make a random walk so 
that after a time t we will have 

((A0) 2 ) r 4((A x ) 2 )A. (49) 

in the case of small Acj> <C 1. Averaging the amplitude of the oscillator over 
the distribution of A0, assuming a Gaussian distribution for A0, the amplitude 
will be reduced by 

where 

1 4 ((A X ) 2 



cos^) ~ e\-^°-f> h) ~ exp ( ) , (50) 

2"r n 



r, 



(51) 



I I -S 



Thus non-specular collisions can be taken into account by the change of the 
damping term in l|32fl : 

l^L + ±. (52) 

Tc T c Tns 

Non-specular wall collisions will thus have a different influence than the gas 
collisions because of the dependence of t w ~ sin a on a. 
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5 Arbitrary Magnetic field geometry 



Our discussion has assumed a magnetic field configuration with G z — dB z /dz 
constant. Pendlebury et al |10| have shown, using a geometric phase argument, 
that regardless of the field geometry the effect only depends on the volume 
average of G z in the high frequency (called by them the adiabatic) limit. In a 
recent note, Harris and Pendlebury have shown that in the case of a field 
produced by a dipole external to the measurement cell, this does not hold in the 
low frequency (diffusion) limit. In this section we discuss this problem using our 
correlation function approach in order to give some physical insight into what 
is happening and display details of the transition from one case to another. 

5.1 Short time (high frequency, adiabatic) limit of the cor- 
relation function 

Reference [TT], has shown that the systematic edm is given, in general, as the 
Fourier transform of a certain correlation function of the time varying field seen 
by the neutrons as they move through the apparatus. Equation (23) of that pa- 
per gives the frequency shift proportional to E as ( oj (t) lies in the x, y plane) 



= -k f d A ^ ■ ^ s ^^(*)>< ^(*-;)] } (53 ) 

W 2 J Q \ +SlliLU T[LU x {t)uj x {t-T)+LU y {t-T)u y (t)} J V > 

It can be shown that the term multiplying sin uj t goes to zero on averaging 
over a uniform velocity distribution ((v x v y ) — 0, v x — v y — 1,2 /%) and using 
V x B — O.Then, for short times, r, 



5u(t) 









COS L0 o T 












COSU> T 







dcu 



LJ (t) X W (t) - — T + - 



1 d 2 U 



duo 



1 d 2 To 



w(<)X (^ r -2dr 2 



(54) 



We are considering values of r so small that the velocity doesn't change in 
that time interval (r < t co «) . 



Then 



duj 
~~dt 
d 2 TS 



(t) = 7 (B xy (t) + !?/c x 
7 ( VB (~x (<)) • ~v 



T-,3 



d 2 B 

dxidxj 



ViVi 
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and 



Sui(t) = - 



2c 



dr cos £ 



- (B xy (t)+ ~vx i?) x f 



9a; 



1 6> 2 I? 



flt T 2 <9t 2 



r 2 + 



(55) 



The term linear in E and r is then 

7 2 " 



Su>(t) = — — / c?r cos o; r 
2c Jo 

= — - / dr cos u) T (ar) 
2 Jo 



V.B • u ) t x ( i? x E 



(56) 



defining 



a = — ^V-B • "w*^ • 1? 



We have now calculated the correlation function for short times. It starts at zero 
at t = and rises as ar. Eventually it will reach a maximum. By concentrating 
on the high frequency (cu ) behavior of Soj the result will be independent of the 
details of the maximum, depending only on a. Thus we can replace ar in l|5()|l 
by sin ar or any function with the same initial slope. Thus we are led to take 



Suit) = — lim / dr 

2 Uo->OcJ Q 



cos uj r ,T sin ar 



lim 



7E 



*co 2 U)„ 



7_E a 



E 

2^m 



VB ■ v 



Introducing components, taking averages and using V • B = this reduces to 



Suigeo — 



-Ev A 



1 



4c£2 \ d 



dB z 



(57) 



in agreement with Eq.. (2) of JI] if, in that equation, R 2 oj^ is replaced by 
(w 2 ) = v 2 /2. We have shown that in the adiabatic (short time) limit the 
systematic (false) edm effect depends only on {^§^-) regardless of the geometry 
of the magnetic field, a result obtained previously by Pendlebury et al ^U] an d 
confirmed in |16| . 

The next order term in (|55|l is easily seen to be of order v 3 t 2 and so will 
average to zero, the next term which contributes will be of order v 4 t 3 and so 
will be negligible in the short time limit we are considering. The condition for 
this to be valid is (vt/L) 2 <§; 1 where L is the scale of variations in the applied 

magnetic field (^3; ~ 
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5.2 Longer time behavior of the correlation function 

For long times the expansion l|54|) is clearly not valid and we must expand in a 
series in the spatial coordinates. We start from 

<W = -I- J dr cos uj t(b' (i)xB' (t-r)^ 

where b — E/c, the brackets represent an ensemble average and 

B' x = B x (t)) - toy 

B' y = B y (V(t))+bv x 



Then we write 
B'(t) x B'(t-T)) =b 



B x ( r (<)) v x (i - r) - u„ (t) B a ( r (f - r)) 
- (B a (V (t - r)) v x (t) - Uy (i - r) fl„ (V (<))) 

(58) 

and expand the held in a Taylor series 



B x (7 (t)) = S s (0,0,0,*) + 



0B a 



l(t) 



0B X 



dy 



d 2 B x 



dx 2 
d 2 B x 



(d 3 B x 



.T 2 (t) + 

y(t)as(t) 
(t) + 



d 2 B. r 



dy 2 
d 2 B 



y 2 (t) 



y(t) 

o 

d 2 B x 



dB, 



dydz 
d 3 B, 



dz 2 



dz 



z 2 (t) 



dzdx 



••(*) 



X (t) z (t) 



dy 3 



y 3 (t) 



(similarly for B y ). Concentrating on the first and last terms in (|58|l and noting 
that there are no correlations between any functions f(xi,Vi) and g(xj,Vj) we 
find 



E 

Xi—x,y 



dB x 



dxi 



Xi (i) 



d 2 B x 



dx 2 



A (t) 



d 3 B x 



dx 3 



4 {t)...)v Xi (t-T)-{(t)^(t-T)} 



where the second term is obtained from the first by interchanging (t) and (t — r) . 

By symmetry we see that (x 2 (t) v Xi (t — r)) = so that the next contribut- 
ing term will be proportional to 



d 3 B x 



dx 3 



(x 3 (t)v x ,(t-r)) 



The first order term will be proportional to 

dB x , dB y 



dx 



dy 



dB z 
dz 
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We see that the condition 



d 3 B. 



R 2 < 



dB. 




Ox 3 



o 



dxi 



o 



will insure that the higher order terms can be neglected. In the extreme case 
considered by Harris and Pendlebury, ^B], this condition is strongly violated so 
our method cannot be applied since the higher order terms remain significant. 

6 Discussion 

We have derived the analytic form of the (2 dimensional) velocity autocorre- 
lation function (vcf) for particles moving in a specularly reflecting cylindrical 
vessel. As the vcf determines the position-velocity correlation function (pvcf) 
we have calculated the Fourier transform of this latter function thus obtaining 
the frequency shift for a single group of trajectories characterized by the angle 
a. Our results duplicate the analytic result found in |l(Jj by direct solution of 
the classical Bloch equations. We then show how the analytic formula for single 
trajectories can be extended to take into account gas collisions and obtain (after 
integrating over trajectories) the linear in E frequency shift for a single velocity 
which has been studied previously, our results agreeing with those obtained by 
numerical simulation of the pvcf with collisions in . We then perform an av- 
erage over the Maxwell velocity distribution using the temperature dependent 
collision times appropriate for He 3 diffusing in superfluid He 4 and obtain the 
temperature dependence of the frequency shift in this case. Figures 4 and 5, 
which are based on an exact average over the Maxwell distribution of the He 3 
velocities, imply that one should be able to control the effect to high degree. 

Due to the heavy mass and slow velocity of the He 3 , Baym and Ebner [T7| 
conclude that the phonon scattering on He 3 is predominantly elastic. Single 
phonon absorption is kinematically forbidden on a single He 3 and can only 
take place as a result of He 3 — He 3 collisions which will be negligible for the 
low He 3 densities considered here. Thus our approach, where we calculate the 
correlation function for an ensemble of trajectories with constant He 3 velocity 
and then average the result over the velocity distribution should be an excellent 
approximation. 

In addition we have discussed the conditions under which the frequency shift 
can be shown to depend only on the volume average of dB z /dz. 



For a general orbit the velocity autocorrelation function is not periodic. How- 
ever, using the form l|15|l for the correlation function, we can calculate its Fourier 



7 Appendix A 



1G 



transform straighforwardly: 



*(a,w) = — A t / cos(u>t)dt + Bi / — cos(ujt)dt . (59) 



(=0,1,2. 



Now ((5 = ojt^) 



(J+1)t u 



cos(wt)<if 



sin((Z + 1)5) - sin(M) 



(60) 



It, 



(i+l)T„ 



i , , , (Z + l)sin((Z + l)5)-Zsin(Z<T) cos((Z + 1)5) - cos(W) 

COS(wi)<K = 1 ; . 

T w OJ LOO 



(61) 

We separate out the terms in (|59fl . according to whether they come from sines 
or cosines in H60(l or (|61|l . 

Taking the sine terms first we have 

((/ + 1) cos(2Za) - Z cos(2 (Z + 1) a)) (sin((Z + 1) 5) - sin(W)) - 

(cos(2 (Z + 1) a) - cos2Za) ((Z + 1) (sin((Z + 1) 5) -lsm(lS)) 

= sin((Z + 1) 5) (cos 2 (Z + 1) a) - sin(Z5) cos(2Za), (62) 

the sine terms in ^ (a, to) are 
,,2 



* s (a, w) = — V (sin((Z + 1)6) cos(2 (Z + 1) a) - sin(Z5) cos(2Za)) 

7TUJ 

/ /(l)-/(0)+ \ 

2 N 2 

= — E(/(' + 1 )-/(0) = — 

■7T7 1 ' * 717 1 1 



1=0 



/(2)-/(l) + 



f(N)-f(N-l) + 
V /(AT+1)-/(AT) / 



2 2 

— / (N+l) = — sin((7V + 1) <5) cos(2 (N + 1) a). (63) 

TXLO 1TLO 



Turning to the cosine terms we have 

2 00 

* c (a, uj) = V S ; (cos((Z + 1)5)- cos(Z<5)) = 



1=0 



2 / N N \ 

>UJ V/=o z=o / 



N 



B N cos{{N + 1)6)-B + ~ B i) cos ( W ) . ( 64 ) 
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where 



B l _ 1 - Bi = 2cos(2Za) - cos(2(7 - l)a) - cos(2(^ + l)a) 
= 2cos(2Za)(l - cos(2a)) = 4 sin 2 a cos(2Za). 



Then 



,,2 



N 



# c (a, uj) = — r 1 - cos 2q + Sjv cos((iV + 1) (5) + 4 sin 2 a cos(M) cos(2Za) , 
^ V U ) 

(65) 

iV AT 

h(S, a) = cos(lS) cos{2la) = - ^(cos(Z (6 + 2a)) + cos(7 (<5 - 2a))) 



2 

i=i /=i 



(66) 



_ 1 / sin ((TV 2a)) sin ((TV + §)(*- 2a)) N 

2 ^ 2sin^ 2sin^l 

except when S ± 2a = 2-7rn in which case /i(<$, a) = iV/2. 
So (we will take the limit as N — > oo) 



sin(Ar+i)a; /siniVx cos ttX \ /sinA^xcos-lx 

v 2 y u_ I 2 i aj„ \ _ 1- / 2 



lim 1 2 = lira = ^— + cos Nx = lim = ^— = 

n^cc sin ^x w^oo y sin =^x J N^-oo \ sin ^x J 

except when x = 2i:m, 

/sin Nxcos ix\ _ , v , x 

lim = 2— = > 2vr(5 (x - 27rm . (67) 

JV^oo I sink / ^ 

x A ' m— — oo 

As a result (|5T)ll becomes (<5 = ujt w ) 



(a, ui) = — ^5 — sin 2 a Y^ (2ttS (ujt w + 2a — 2-Km) + 2irS (ojt w — 2a — 27rra)) . 

TTU) Z T W L — ' 

m — — oo 

(68) 

Using S (ujt w + 2a - 2vrm) = ±5 (w + 2a ^ m ) we have finally 



2w 2 sin 2 a /_ / 2a — 27rm\ . / 2a + 2nm 



m= — oo 
9-2 oo 

'i^ cin 



_ , , 4V sin a v— * I I ^a — civiii \ . / 



2 "S^ E («(«-«+(«)) +%-^(a))), (69) 



m— — oo 



where w=(a) = v r 
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8 Appendix B 

Using the formula [TH| 



cot 6 



I H ! ^ \ 1 




oo 1 



n— — oo 



n(6 — rar) 
(^2 means n = 0, excluded) one gets: 

cot 6? — cot <p 1 

os-61 ~ ^ (7m-<zj)(7rn-0)' ^ ^ 

r n— — oo v ' x 7 

After differentiation with respect to tp we obtain: 

\p 1 = 1 cot 9 - cot y? 

„^(™^) 2 (™-#) (y-^sin 2 ^ (v>-0) 2 ' 1 j 
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9 Figure Captions 

Figure 1. Note the curves are for a single fixed velocity. The velocity dependence 
is contained in the normalization of the frequency scale, ui r = v / R. 
Figure 2. Trajectory of a particle in a cylindrical cell. 

Figure 3. Normalized frequency shift for a constant velocity as a function of 
normalized applied frequency, to' = lo R/v, for different values of the damping 
parameter r Q = R/X. Solid curves - results of the analytic function, equations 
<|43l) and (|44|l . Dotted lines numerical simulations from ref. [TT]. red - r D = .2, 
green - r D = .5, cyan - r Q = 1, black r G — 2, blue - r Q = 4, magenta - r D = 10. 
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Figure 4. Normalized velocity averaged frequency shift vs. reduced fre- 
quency u>*=u R/ j3 (T) for various temperatures using the temperature-dependent 
mean free path for He 3 in He 4 . 

Figure 5. Normalized velocity-averaged frequency shift, ^ (uj*,T) vs. tem- 
perature, T, for various reduced frequencies u>*=u> R/ f3 (T) using the temperature- 
dependent mean free path for He 3 in He 4 
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